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ABSTRACT 

We show that radiative feedback due to reionization has a pronounced effect on the extent of 
mechanical feedback due to galactic outflows. The photoionization of the Intergalactic Medium (IGM) 
suppresses low-mass galaxy formation by photohcating the gas and limiting atomic line cooling. The 
number of low-mass galaxies is central for the enrichment of the IGM as these objects have the capacity 
to enrich a significant fraction (by volume) of the Universe. We use a modified version of our galactic 
outflow model, combined with a simple criterion for suppression, to investigate the potential impact 
upon the IGM. We find that this suppression strongly reduces the enrichment of the IGM and is 
sensitive to the reionization history. We also investigate the contribution of halos of different masses 
with varying degrees of suppression. 

Subject headings: cosmology: theory — galaxies: formation — intergalactic medium — galaxies: dwarf 



1. INTRODUCTION 

Feedback is an essential component of our understand- 
ing and modeling of the intergalactic medium and galaxy 
formation. This feedback falls broadly into two cat- 
egories: radiative feedback and mechanical feedback. 
They both act to deposit energy into the intergalactic 
medium (IGM) which in turn suppresses star formation. 

Radiative feedback is dominated by ionizing photons 
and their capacity to reionize the Universe, keeping it 
ionized until the present epoch. This ionization pho- 
toheats the IGM and suppresses a crucial source of 
cooling: atomic line cooling. Radiative feedback acts 
to suppres s the formation of low-mass galaxies (e.g . 
Reel fl98l lEfstathioul H99l iThoul fe Weinberg! [1996 : 



Navarro, fc Steinmetd 119971 : IKitavama. fc Ikeuchil boOO : 
Dijkstra et all 12004 hereafter DHRW) since the gas is 
unable to cool and collapse into the dark matter halos to 
form stars. 

Mechanical feedback results from a variety of sources, 
but up to z k, 2 the main source of this feedback is 
thought to be galactic outflows caused by many co- 
eval Type II supernovae exploding during the starburst 
phase. Population III stars may also contribute signifi- 
cantly but early result indicate they will be sub-dominant 
(|Norman. O'Shea, fc Paschodl2004) . These outflows do 
not simply deposit energy into the IGM but also deposit 
metals. These metals increase the gas cooling efficiency 
and also pr ovide a convenient observati o nal tr acer of this 
feedback. IScannapieco fc Broadhurstl ((2001) nn d that 
most mechanical feedback is caused by low-mass galax- 
ies, however, this is without the inclusion of radiative 
feedback. 

In this work we seek to investigate the degree to which 
radiative feedback can suppress the formation of low- 
mass galaxies that form in dark matter halos with virial 
temperatures > 10 4 K (below which atomic cooling is in- 
efficient). Hence we do not consider the contribution to 
our outflows of minihalos and stars which form with a H2 
cooling mechanism. We do not make specific predictions 
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for the degree of suppression and so the extent of metal 
enrichment, rather we constrain the potential impact of 
this effect. We use analytic models for the expansion 
of gala ctic outflows described in lPieri. Martel. fc Grenonl 
(2007) (hereafter PMG07), which self-consistently in- 
clude the impact of ram-pressure stripping of halos and 
metal enrichment. We make a simple parametrisation of 
the degree of photoionization suppression by choosing a 
step-like reionization and using a suppression criterion 
taken from DHRW. 

This paper is set out as follows. In [J2]we describe our 
criterion for photoionization suppression. In [J3]we sum- 
marize the model of galactic outflows used. The results 
from our simulations are set out in SJ3J We discuss the 
implications of our results and conclude in iJS] 



2. THE PHOTOIONIZATION SUPPRESSION CRITERION 

Our criterion for photoionization suppression can be 
expressed as follows: halos ionized before turnaround, 
with a mass below a certain threshold, fail to cool a 
'sufficient' quantity of gas in order to form a galaxy. 
This mass threshold corresponds to a particular halo cir- 
cular velocity v c \ IC and a particular virial temperature 
?vir- We specify this quantity of gas in two ways: 1) 
Wcirc = W1/2, where half the halo baryons or more cool, 
and 2) u C i rc = vq, where galaxy formation is only sup- 
pressed in halos where no gas cools. The former is an 
acceptable approximation for our Monte Carlo simula- 
tions and is used in our fiducial model. The latter is also 
considered in the interest of placing a lower limit on the 
suppression. 

In Figure [T] we use the results from DHRW's Figure 3 
to derive a fit to the evolution of the critical circular 
velocity for a given collapse redshift, z co ii- We show 
their results for three cases: t>x/2 f° r an evolving UV 
background amplitude (J = J21 x 10 -21 ) taken from 
IKitavama. fc Ikeuchil ((2000), V1/2 for a fixed amplitude 
J21 = 10~ 2 , and wo f° r an evolving amplitude. All 
three have been fitted with the straight lines shown us- 
ing the least-squares method. Where an evolving UV 
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Fig. 1. — The threshold halo circular velocities prevented from 
forming stars for halo collapse redshift. All points are taken from 
DHRW. The squares and the crosses show the case where v c - lIC = 
Di/2- The squares are for an evolving UV background and the 
dashed line is a least-squares fit to these points. The crosses show 
the case were the UV background amplitude has a fixed value of 
J21 = 10 — 2 and the dotted line is their fit. The dot-dashed line 
is the fit to the case where v c i IC = vq and the UV background 
evolves, given by the triangles. 

background is used we find circular velocities, 

t>i /2 = [89.4 ±2.4- (6.3±0.3)z coU ]kms- 1 , (1) 

which corresponds to our fiducial model. We also con- 
sider modifications to this model using 

v = [42.1 ± 2.5 - (2.6 ± 0.4)z coU ]kms" 1 (2) 

and for a UV background with a fixed amplitude, 

v 1/2 = [69.2 ± 2.4 - (4.3 ± 0.3)z con ]km S - 1 . (3) 

These dependences are used in the following work 
to provide a criterion for photoionization suppression. 
These circular velocities can be related monotonically in 
a straightforward manner to a virial temperature, mak- 
ing the criterion simple to incorporate into our cooling 
time calculation (note that for T vu - = 10 4 K, u C i rc w 
17km s _1 ). We require that the virial temperature is be- 
low the value derived for suppression to occur. DHRW 
argue that the choice of evolving J21 constitutes an upper 
limit on the expected v c irc but it is notable that the ef- 
fect on mechanical feedback of switching from equation 
(Q} to equation ([3]) are almost indistinguishable in the 
results that follow. 

DHRW perform their simulations using a one- 
dimensional spherically symmetric hydrodynami cs code 
based on that used by iThoul & Weinberg! (|1995l ). They 
consider a halo which is crossed by an ionization front 
during collapse; this occurs at z = 17 in their simula- 
tions. We prefer to choose the crucial reionization red- 
shift, which we call z r . We do this in our fiducial model 
by assuming that, so long as the halo itself is ionized be- 
fore its turnaround redshift, 2; tU rn, the precise redshift of 
reionization is unimportant. This is a reasonable approx- 
imation but will underestimate the degree of suppression 
at high redshift since equations (HI)-© are derived from 
DHRW and for their z co n = 11, z tuin = 18 and at these 
redshifts z r < z tU rii- As result we also show the other 
extreme case, allowing suppression when the z r > z co \\. 

It should be noted that our choice of z r is not intended 
to be characteristic of reionization as a whole, but rather 
a good representation of when low-mass halos are ion- 
ized. We do not select one value for z r but rather use 
upper and lower values of z T — 11 and z T = 6. The 



higher redshift value is chosen to correspond to the red- 
shift of reionization derived from the 3-yr WMAP data 
relea se in a step-like reionization scenario ()Page et al.l 
2006). z r = 6 constitutes a conservative lower limit on 
the significance of photoionization suppression based on 
the end of the reionization epoch which may have been 
detec ted in quasar absorption spectra (e.g. iFan et ahl 
l2006HBecker. Rauch. fc Sargent II200I . 

3. THE MODEL 

For a complete description of the model used see 
PMG07. We will describe some of the main features 
of the code and those particularly significant for this 
work. The approach is a Monte Carlo simulation of a 
(12/i -1 Mpc) 3 box, which constitutes a cosmological vol- 
ume down to the end point of our simulations: z = 2. A 
ACDM cosmology consistent with ISpergel et aT~l (2006) 
is used and a Gaussian random realization of struc- 
ture is formed on a 512 3 grid. The dark matter ha- 
los, their mass, M , and their collapse redshifts, z C oii; are 
found. Mergers are included although only in a non- 
dynamical manner. Upon collapse the gas inside these 
halo s virializes and begins to co ol. We use MAPPINGS 
III 2 ([Sutherland fc Dopital [T993) to calculate the cooling 
time. 

When cooling is complete galaxies are formed with a 
star formation efficiency /». Type II supernovae then 
inject energy into the IGM, preferentially taking a path 
of least resistance out of large-scale structures with an 
opening angle, a. In our fiducial model we use isotropic 
outflows but also consider an anisotropic case. The ex- 
pansion of these outflows out of the gravitational poten- 
tial well of the host halos is calculated using the thin 
shell approximation. We assume that the IGM is ion- 
ized throughout the expansion of the outflows. This is 
an acceptable approximation since the epoch of outflow 
expansion is later than most (if not all) of reionization 
and so (pre-turnaround) halo suppression. 

4. RESULTS 

We now have the modifications required to perform our 
simulations of mechanical feedback with a simple descrip- 
tion of photoionization suppression of low-mass galaxy 
formation. We run our simulations to their end point of 
z = 2 and compare the regions which have experienced 
galactic outflows with the underlying gas density field. 

As in PMG07 we describe this IGM density field, p/p, 
by producing a realization of gas structure using Jeans 
filtering. We map the resultant gas density field to a 
lognormal probability distributio n function to reproduc e 
a degree of non- linear behavior (|Bi fc Davidsenl fl997f) . 
Every point in the N = 512 3 grid is compared with the 
outflow locations to determine if it is within the volume 
of an outflow, the gas density of that point is also noted. 
In Figure[2]we see the results of this statistic. Each panel 
shows the cumulative statistic (N' pl<p /N) of the number 
of grid points enriched at a given density N' , which gives 
the volume filling factor of enriched systems below this 
density. This is shown for the case were no suppression 
of gas infall is included and the cases where it begins at 
z r = 11 and 6. The second and third panels show hard 

2 http:/ /www. mso.anu.edu.au/~ralph/map. html 
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Fig. 2. — The volume filling factor of enriched regions with a gas overdensity < p/p. All panels show the statistic when no photoionization 
suppression is used in our model with isotropic galactic outflows (solid line). The first panel shows the fiducial model of photoionization 
suppression described in the text for z r = 11 (dashed line) and z r = 6 (dotted line)). The second and third panels show the cases where the 
fiducial model is modified to suppression when z T > z co n and v c - lTC = vo respectively. The fourth panel shows the case where our fiducial 
model is modified to include anisotropic outflows with a = 100°. Here the dot-dashed line shows the statistic with anisotropic outflows and 
no suppression. 

fcirc — V Q an< i Zr = 11 are used and the volume enriched 
rises for halos of mass ~ 1O 9 5 M . This is mainly be- 
cause the lower circular velocity threshold allows such 
halos to escape photoionization suppression. In the sec- 
ond case, anisotropic outflows are used with no suppres- 
sion resulting in a peak at around 10 9 M Q . This results 
from a combination of factors related to clustering but 
mainly: a) a reduction in stripping since outflows take a 
path of least resistance away from other halos, and b) a 
lower volume per outflow combined with less overlap in 
outflows. 

It is striking that suppression of gas infall in the z r = 6 
case makes little change to the impact of halos of mass 
< 10 8 5 Mq in all but the strongest suppression (shown in 
the second panel). This is because, in a ACDM model, 
most halos in this mass range turn around early, before 
z T . This also explains why the third panel of Figure [2] 
shows only a small change in the z T = 6 case compared 
to the first panel: the significant lowest mass halos are 
unaffected and so the results are insensitive to the v c i TC 
threshold used. The z r = 11 case shows a significant sup- 



upper and lower limits (respectively) on our suppression 
criterion and its impact. 

The first panel shows the fiducial case with a suppres- 
sion criterion using equation |T]) , ztum < z r and isotropic 
galactic outflows. This shows that photoionization feed- 
back has a pronounced effect at all densities and may re- 
sult in a ~ 75% reduction in the volume of the Universe 
affected by galactic outflows. Even in the conservative 
case of z T = 6, a ~ 25% smaller volume is enriched. This 
result uses the condition that a galaxy can only be sup- 
pressed when its host halo has turned around after it has 
experienced reionization. If we require instead that sup- 
pression may occur even after turnaround and up until 
collapse, we obtain the results shown in the second panel. 
This certainly results in excessive suppression. 

The third panel shows the case where we use equa- 
tion |2j) to provide the circular velocity for which sup- 
pression occurs. This is rather conservative as suppres- 
sion only occurs where no gas falls into the halo; all other 
halos arc unaffected. The case where z r = 6 is almost 
unaffected by this adjustment, while the degree of sup- 
pression in the z T = 11 case is reduced from ~ 75% to 
~ 60%. 

We show an example of the impact of this suppression 
for anisotropic outflows in the last panel of Figure[2l This 
includes a run with an opening angle, a — 100° for no 
suppression, z x — 6 and z T = 11. Also shown is the case 
of isotropic outflows and no suppression for comparison. 
This illustrates the difference in the nature of the two 
effects: anisotropic outflows have a strong impact on the 
degree of enrichment in low density regions (see PMG07) 
while this suppression affects all densities equally. 

Figure [3] shows the breakdown of volume filling factor 
enriched by different halo masses for all the models shown 
in Figure [21 Note that the sum of these histograms is not 
the total volume filling factor as some regions have been 
enriched more than once by galaxies of different mass. 
In all cases the dominant sources of mechanical feedback 
are galaxies in halos of mass 10 8 M Q < M < 10 9 M o . It 
is clear that low mass halos dominate the volume filling 
factor of metals, even when photoionization suppression 
is included. 

In all but two of the cases considered, the volume en- 
riched peaks around 1O 8 5 M and decreases monotoni- 
cally with increasing mass. In the first of these two cases, 



pression of galaxies in halos of mass < 10 Mq, however, 
the impact of suppression on these halos is still small 
compared to the impact on higher mass halos and they 
continue to dominate the overall volume filling factor. 

5. DISCUSSION AND CONCLUSIONS 

Our analysis of the influence of photoionization sup- 
pression of galaxy formation is based on a simple pre- 
scription taken from the findings of DHRW. We use an 
analytic description of galactic outflows in a Monte Carlo 
cosmological simulation (PMG07) . In our fiducial model 
we have chosen to suppress the formation of galaxies in 
halos which allow the infall of half their baryons or less. 
This occurs only when the halo has experienced reion- 
ization before its turnaround redshift. These choices, as 
well as the use of a step-like reionization, are justified 
simplifications, but in each case we have also bounded 
the simplification with the lower-limiting case. In all our 
result we find that this form of radiative feedback signif- 
icantly reduces the physical extent of enrichment. 

The largest unknown in our prescription for suppres- 
sion is the reionization redshift. Work must be done to 
determine what redshift satisfactorily characterizes when 
halos which dominate the enrichment of the IGM are 
reionized. This is further complicated interdependence 
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Fig. 3. — The volume filling factor (Nij/N) of regions which have been enriched by galaxies in halos of mass, M. The lines correspond 
to the same models as shown in Figure [2] 



of rei onization history and photoionization suppression 
(e.g. Illiev et all 120071 ) . Also the values taken for v c irc 
represent upper limits since self-shielding is not included 
(see DHRW) . It may be necessary produce models which 
describe both radiative transfer and galactic outflows in 
order to address this issue with any precision. 

We have shown that this suppression has a pro- 
nounced impact on the extent of galactic outflows, 
but one should note that changes in the star for- 
mati on efficiency result in a similarly pro nounced ef- 
fect (jScannapieco, Ferrara. fc Madaul 12003 ) . We chose 
a value of /» = 0.1 which is widely accepted (e.g. 
iBarkana fc Lo eb 2000) , but we find that doubling (halv- 
ing) /* leads to ~ 50% larger ( ~ 67% smaller) volume 
filling factor in all our simulations. This change in /» 
would of course modify the reionization history (which 
would counter the impact of changes in _/*) but we have 
not considered this here. 

It should be noted that the non-dynamical nature of 
our models and the way they treat mergers means that 
we will tend to have more low-mass halos than one would 
expect. As Figure [3] shows, one may not simply argue 
that fewer low-mass halos would result in a weaker sup- 
pression signal. The least massive halos which have a 
significant impact on the IGM by z — 2 are less affected 
by this suppression as they are more likely to turn around 
before reionization. 

Observationally the volume filling factor of the en- 
riched IGM is rather poorly constrain ed but it appears 
that the Universe is > 5% enriched (jPieri fc Haehneltl 
2004). If we assume that the star formation efficiency is 
not significantly higher than our chosen value, our results 
seem to rule out the strongest photoionization suppres- 
sion in our mo dels. Metals have b een observed in the 
IGM at z > 6 ([Becker et al.l [2006T) and our results are 
further indication that metal enrichment and reioniza- 
tion are no-longer two distinct fields but part of a grow- 



ing field of reionization chemodynamics. 

Using our fiducial model, we find that photoionization 
suppression may result in a ~ 75% reduction of the frac- 
tion of the Universe enriched at all densities in the IGM. 
Furthermore, we find that even using a variety of con- 
servative modifications to the suppression, we still find a 
~ 25% reduction in the enrichment of the IGM. We also 
find that the least massive halos (< 10 8 5 Mq) are almost 
unaffected by late reionization and even with our best 
reionization redshift we find that these halos are sup- 
pressed to a relatively low degree. Despite this suppres- 
sion, low-mass galaxies continue to dominate the metal 
enrichment of the Universe by volume. 

Further work is required to determine the precise de- 
gree of suppression and the largest unknown in our mod- 
els is the local reionization history experienced by halos. 
Also a focused study constraining the amount of star for- 
mation in halos of various masses which are irradiated 
at varying stages of collapse would significantly improve 
our analytic prescription of this effect. Observationally, 
30m class telescopes are required to significantly improve 
constraints on the volume filling factor of metals. In the 
meantime, more studies must be performed investigating 
the nature and environment of enriched systems, while 
extending QSO absorption studies to higher redshift in 
order to study the evolution of metallicity at the end of 
the reionization epoch. 
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